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Abstract 

We prove the stability of the one dimensional kink solution of the Cahn-Hilliard 
equation under d-dimensional perturbations for d > 3. We also establish a novel 
scaling behavior of the large time asymptotics of the solution. The leading asymp- 
totics of the solution is characterized by a length scale proportional to fs instead 
of the usual 1 2 scaling typical to parabolic problems. 

1 Introduction 

The Cahn-Hilliard equation is a fourth order nonlinear evolution equation for a real 
valued function u{x, t) defined on some spatial domain x G il C M'^: 

dtu = A{-Au-lu+^u^) (1) 
u(x,0) =.g(x), xeO. (2) 

The nonlinear term inside the brackets in the RHS has three zeros, u ^ 0, ±1 where 
the first is linearly unstable whereas the others are linearly stable. 

The CH equation is used to model phase separation in mixtures of two substances 
A,B (binary alloys) so that u(x, t) describes relative concentration of the substances 
and the zeros ±1 correspond to pure phases of A or B. 

When a random initial condition g is given the solutions of Q typically exhibit in 
numerical simulations phase segregation, i.e. domains of phase A and B start to form 
and increase in size until they reach sizes comparable to the domain size. To understand 
such extensive behavior of the solutions it is natural to consider Q in the whole space 
= M'' which will be assumed in the present paper. 

In one dimension a single phase boundary is described by a stationary solution of 
([0, the so called kink solution. This remains a solution also in dimensions d > 1 and 
is given by 

uq{x) — tanh(ix) (3) 

where x is the first coordinate of x. Thus up to exponentially decaying tails, Q de- 
scribes a situation where we have phase B in the domain x < and phase A in the 
domain 2: > 0. 

The presence of the fourth order derivative in Q makes the mathematical analysis 
of the CH equation much harder than analogous second order equations. The absence 
of a spectral gap in unbounded domains due to the Laplacian multiplying the RHS 
also complicates matters. In |2| and 13 the stability of the kink solution in one dimen- 
sions was proved. Moreover in ||2|| the following leading asymptotics for u{x,t) was 
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established: 



Ad/,. _xi\ B d _xi , 

u(x + a,t) = uo{x) ^ p— uo(x)e « H p-;— e « + ) 

Vt d.T V / da; 

(in sup norm) where the constants a,A,B depend on the initial data and the function 



^uo(a;) = (2cosh2(ix))-i (4) 



decays as 2e 1^' . Thus, for large times one observes a translated front (from the origin 
to a), a perturbation of size 0{-^) localized near the origin and a perturbation of size 

0{j) extending to an interval of size ^/t around the origin. The latter exhibits typical 
diffusive scaling between space and time. 

In the present paper we prove stability of the kink solution when the spatial di- 
mension d > 3 and establish the following asymptotics for it. Let us agree to denote 
variables in M'' by the letters x or y with x = (x, x) € K x M*^"^ and so on. Moreover, 
k or p will only be in M'^^^ with k := |k| and the same for p. Define the functions 

and for t > 

ci,{^,t) = t-'^r{^) (6) 

or in terms of Fourier transform, 

0(k,t) = 0*(t^k) ^e-^*!"^!'. (7) 
Let the initial datum (|2j be given by 

g{x) = uo{x) + h{x) (8) 

and the function h satisfy 

\\h\\x ■.= sup\h(x)\{l + \x\Y <5 (9) 

X 

where r > d+ 1. Then we prove 

Theorem 1 . Let d > 3. For S small enough the equation Q has a unique classical 
solution satisfying for t > \ 

u{x, t) = uq{x) + —dxUQ{x)(t>{y., t) + u(x, t) (10) 

where 

swp\u{x,t)\<Ct~^t~^ (11) 

X 

and A ~ J^^ h{x) dx. 
Remark 1 . Since 

A A d-i 

Mx) + -a^uo(a;)0(x, t) = uo{x + -0(x,t)) + 0{t~^ — ) 
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we see that (I10> describes a a front that is translated in a domain of size around 

_ d-l 

the origin by a value of the order t s . Note that in contrast to one dimension the 
translation of the front tends to zero as time tends to infinity. This is because a localized 
perturbation is not able to produce a constant shift in the whole transverse space M"^^^. 
However, the perturbation does not decay in the standard diffusive fashion but with the 
different power of time: >/t is replaced by i 3 . This scaling was argued to be present in 
the linearized CH equation in (S) and Q. We prove actually more detailed properties 
on the spatial behavior of u, see Proposition^ 

Remark 2. In two dimensions the nonlinearity becomes "marginal" in the terminology 
of (Q. We do not know whether the asymptotics proved in the Theorem persists there. 

The remainder of the paper is organized as follows. In Section|2]we present how the 
problem is reduced to a nonlinear parabolic Cauchy problem with small initial data. We 
also state the main estimates for its semigroup kernel needed for the nonlinear analysis. 
In Section|3lwe use these estimates to bound the nonlinear terms in the integral equation 
corresponding to the above mentioned Cauchy problem, thus proving the main result. 
The proofs for the crucial semigroup estimates are presented in Sections l4lj 101 

2 Linearization 

We start by separating the kink solution 

u = uq + V. (12) 
Recalling that uq solves the CH equation we get for v the equation 

dtv =^ Cv + A{^uav^ + ^v^), (13) 

where 

C~~A{A+^-lul). (14) 

This linear operator will play an important role in the analysis since it will pro- 
vide the leading asymptotics. Indeed, we will solve the equation (I13> with the initial 
condition vq ^ hhy studying the equivalent integral equation: 

v{t) = e*^h+ f e^'-'^^ A{^uov{sf + ^v{sf) ds. (15) 
Jo 

Let us therefore discuss the properties of the semigroup exp{tC) generated by C. 
Write the operator C as 

C = AH with H:=-A+l + V and V{x) := cosh{^)-^ (16) 

Since C is constant coefficient in the transverse x direction it will be convenient to 
work in a mixed representation, with Fourier transform in these variables. Thus given 
f : M.^ ^ C denote by f{x, k) the Fourier transform with respect to the d — 1 last 
coordinates. In this representation —C becomes 

{-AHfrix,k) = iDkHkf){x,k) (17) 

where 

Dk = -dl + e, Hk = Dk + l + V 
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Figure 1: The spectrum of DkHk (for small k). 

and we denoted |k| by k. From now on we will work in the (x, k) representation and 
for notational simplicity drop the hats from Fourier transforms. 
The semigroup is then written as 

(e*^/) (x, k) = / dy K{x, y, k, t)f{y, k) (18) 
Js. 

with K{x, y, k, t) the integral kernel of the semigroup of the operator —DkHk- In this 
notation the integral equation il5\ becomes 

v{x, k, t) = (e*^/i)(x, k) + ^ J^^y^ k^)K{x, y,k,t- s) 

■ (|Mo(y)w*^(2/, k, s) + k, s)) ds dy 

(19) 

where the Laplacian was integrated by parts to act on the semigroup kernel and * 
denotes convolution in the k variable. 

We will express the semigroup as a Dunford-Cauchy integral of the resolvent ker- 
nel: 

K{x,y,k,t)^ f p,e-'^\C-DkHk)-\x,y) (20) 

where F is a suitable curve around the spectrum of —DkHk- 

The resolvent kernel in i20\ may be studied by standard ODE methods as was done 
in the one dimensional case in |2|. This is rather straightforward but tedious and in 
this section we motivate and present a lemma summarizing the estimates needed for 
the nonlinear analysis. The proof of the lemma is given in Sectionsl4ljlOI 

The spectrum of the operator DkHk is on the positive real axis. Furthermore, there 
exists a fco > such that for k small, k < ko, the spectrum contains an isolated point 

Co-ifc3 + 0(fc4) (21) 

and the rest of the spectrum is on the semiaxis § = [|fc^,oo), see Figure [J The 
resolvent has a simple pole at Co and is analytic in the complement of {^o} U §. For k 
larger than ko the spectrum lies in [6, oo) with 6 > 0. 

Since the function V in ( I16> decays exponentially (as — 6e^l^l) for large x the 
behavior of the resolvent for x and y large is determined by the functions in the kernel 
of the constant coefficient operator 

C-Dk{Dk + l) 

obtained by setting V to zero. These are given by e^^' with C — (— /i^ + k'^){—ii^ + 
A:2 + l) = 0i.e. 

M = ±v/i + P±i^i + 4C. 
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Figure 2: Integration paths for small values of k. 



For large times t the main contribution to ( I20> comes from small C,. In that domain the 
eigenvalues are approximately 

Ml « Vl + fc' + C (22) 

and 

~ Vk^-C (23) 

and their negatives. 

The integration contour in ( I20> will be chosen as follows. Let first k < kg. Then 

(a) For k < the contour is as in Figure |2(a)| In the neighborhood of the origin 
the eigenvalues are /ii « 1 and Rc /i2 ~ ct^^ i.e. the decay of the resolvent is a 
combination of 

g-(i+o(t-^))|x| g-o(t-4)|.| 

(b) For < k < ko the contour circles the pole and the semiaxis § as in Figure |2(b)| 
At the pole the eigenvalues are 

fil=l + 0{k^), I^L2^ {l+0{k^))k. 

On the second part of the contour and \(\ < Ck^ fii is close to 1 and fi2 has real part 
0{k) i.e. larger than t^^. 

The resolvent has a representation in terms of the functions in the kernel of C and 
its adjoint. At fc = ^ = these are explicit and for small k and C they may be studied 
perturbatively. We will need explicitly a few of the leading contributions to K for small 
k. The following lemma summarizes this. 

Lemma 1 . There exists fcp > such that for all k < kg, t > 1 and all x, y the integral 
kernel K{x, y, k, t) of the semigroup of the operator ^D^Hk may be decomposed as 

K = Ko + Ki (24) 

where 

Ko{x,y,k,t) = ^duo{x)Z{y,k,t) (25) 
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and 



Ki{x,y,k,t) ^ ^=(e -e ' \l^y^oe + K2{x,y, k,t). (26) 



The function Z is even in y and has the property 

\ f Z{y,k,t)h{y,k)dy-e-i"'" f h(y)dy\ < C\\h\\xt-h-i"'" (27) 

whereas 

\K2{x,y,k,t)\ 

< fct-e--l^-l ifk^t-K 

We will decompose in the same way also the kernel occurring in eq. il9\ : 

S := [dl - k^)K = So + S^. (29) 
The following estimates hold: 
Lemma 2. (a) Let k < t^^/^ andt> 1. Then 

\So{x, y, k,t)\< C(t-5e-5(l-l+l^'l) + t-ig-^d-l+Ti)) (30) 

and 

\Si{x, y, k,t)\< Cit-^e-i^'^-y^ + t-ie"-*-'''!-!--!!'! + t-^/^er'^). (31) 
(b) Let <k<ko andt > 1. Then 

\Soix,y,k,t)\ < C(fce-^(l^l+l^l) + fc2e-5l^l-^'^-l2'l)e-^*'''' 

+ (^(i^se^^d^l+l!/!) + i-ig-5l^l-cfc|!/l)e"^*'=^ 

and 

\Si{x,y,k,t)\ < C(fc^e-"''l^"^l +A:3e-"'=I^I-"I^I +fc2g-i|:.-y|)g-itfc^ 



- 2 ifc^ 

(33) 

Remark. In order to get a feeling for the various terms, the following intuition is 
useful. For large t and x and for small k the resolvent is built out of functions that 

decay approximately as e"'^', e"*^'^' and e~* ^ '^L Each x derivative brings a factor 
of 1, k and t~2 respectively. Thus, e.g. the third term in ( 13 U behaves as the second 
derivate of e* ^ . 

Finally, the large k or short time behavior of the semigroup is dominated by the 
fourth derivatives in the symbol; 

Lemma 3. (a) Let k > ko and t > 1. Then 

\K{x,y,k,t)\, |5(.T,2;,fc,t)| <Ce-5'='*e-^l^-^l (34) 
(b) Let t < I. Then for all k 

\Kix,y,k,t)\ < Ce-^^'H--^e-*~^\''-y\. (35) 

and 

\S{x,y,k,t)\ < Ce-^^'H-ie-^~^\''-y\. (36) 



6 



3 Proof of the Theorem 



We solve ^ using the contraction mapping principle in a suitable Banach space. For 
each t > 1 we define the Banach space Xt of continuous functions / : M x M'^^^ C 
as follows. First, let m = r — + 1 (r is defined in (|9} so m > 2) and define 

For t > 1 let 

kt = min{fc, 1} + ^ (37) 

and the same notation is used for any positive real number in place of k. 
The norm in Xt is defined to be 

WfWt := sup{iuix) + kMhx)y\l + khr\f{x,k,t)\ (38) 

k>0 

Here n can be taken arbitrary number larger than The main estimate is the 

following 

Proposition 1. There exists a 5 > Q such that if the initial data satisfies |/i(x)| < 
5{\ + |x|)^'' then the equation has a unique classical solution v such that, for all 

t >1, v[t) £ Xt and 

\Ht) - ^a,uo(a:)0(x,O||t < C5r^. (39) 
Remark. Since for ijj & Xt 

|V^(x)| < j Ak{uj{x) + ktUj{ktx)){l + kh)-''\\^\\t<Ct-^-^\\^\\t 

the sup norm of v{t) is bounded by Ct 3 A and the Theorem follows. 

The rest of this section is devoted to the proof of Proposition ^ The proof splits 
into short times and long times. For short times we have the following lemma (see (|9} 
to recall the definition of || • ||x): 

Lemma 4. There exists a 5 > Q such that if the initial data satisfies |/i(x)| < (5(1 + 
|x|)^'' thenv{l) e Xi, 

\Hl)\\x<C5 and \\v{l)\\i<C5. (40) 

Proof. This is quite standard: the leading symbol of the linearized equation is smooth- 
ing and preserves polynomial decay. To prove the second estimate we need to control 
the large k behavior, in view of the definition of the norm in ( I38> . Here it is more nat- 
ural to work in the x representation and derive sufficient estimates for the derivatives. 
Hence let X^p' be the space of p times continuously differentiable functions ^ C 
with the norm 

ll/llx(p) — max||a^/||x 
ImI<p 

'The third power comes from 1271 . The limit will provide sufficient k-integrability for the proof of 
Lemmal^ 
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where /i is a multi-index. We proceed to show that if the initial condition is in X^p") 
after an arbitrarily short time the solution will be in X^p~^^\ 
Write (I13> in the form 

dtv:=^-A^v + AN{v) (41) 
with the initial condition v{0) = h and 

N{v) := i((3Mo - l)v + Suav^ + v^). 
This is equivalent to the integral equation (after integration by parts) 

v{t)^er*^'h+ / (Ae-(*-")^')7V(w(s))ds, (42) 
Jo 

which can be differentiated: 

pt 

d^v{t) = e~*-^'d^h+ Ae^(*-'')^'a^iV(w(s))ds. (43) 



From the explicit Fourier integral representation it is easy to see that the integral kernel 

G(x, y, t) of e~* ^ satisfies the bound 

|9;G(x,y,t)| < C,t-^e-''"'"\^-y\ 
for all multi-indices i/ and x, y e M.'^. Thus for any f G X 

ll/llx /■ 1 , , cii/iu 

<GSUpe 2*1/4 _ j— J— / —7776 2tl/-l dy < 



(i + |y|)'-7E<<i^/4 ^-(i + |x|)'- 

when t < 1 and for any g G C([0, 1], X) 

*Ae-'*-)^^g(.)d.< ^ sup |lg(.)||^ f \ ds 

(1 + |x|)'- ,,g[o.i] 

It follows that for any small enough t > ( I42> can be solved by the contraction map- 
ping principle in the Banach space C([0, t], X'^p^) with the maximum norm. Namely, 
there exists Sq > s.t. if ||ft.||_Y(p) < ^o then ||i'||c([o,t],x<p)) — C'll^llx'?)- 
Differentiating ( I43> yet again yields 



Estimate the integrals as in the previous case to get 

|a,^aMt,x)| < (1 ^^^1). (^ll^x^^lU + sup \Kv\\cao,r],x) 

Thus the solution gains another derivative in an arbitrarily short time. This can be 
iterated to get the estimate |K'(1)ILy(3") — C'II^ILy, which proves the lemma. □ 
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Let us denote 

f = vil) 

and write our integral equation (I19t for t > 1 

+ f f S{x,yXt^s){luo{y)v*\y,k,s) + ^v*''{y,k,s))dsdy (44) 



We want to prove that v{t) satisfies the estimate i39\ . Let us start with the solution of 
the linear equation: 

Lemma 5. Let t > I. The solution to the linearized problem is given by 

e^'-'^^f = vo + vi (45) 

vo{x,-k,t)^^duo{x)e-i"'' (46) 

with A = Jjj^ / and 

IK'i|lt<Ct-^(|l/|U + 11/110. (47) 
Proof. Let first t > 2 and k < ko. Using (I24> decompose 

By (|25} and ^ 

A 13 

wo(i,a^,k) = —duo{x)e~s'' * + w{x,'k,t) (48) 

with 

/ /(x)dx 
Jr'' 

and 

|i&(:E,k,t)|<C|l/|lxt-^e-3^-'*e-l-l (49) 
(we used |i9uo(x)| < 2e~l^l) whereby 

ll^&lk <Ci-3|l/|lx. 
For wi we use the decomposition of eq. ( l26l to write 

Wl = Wii + Wi2. 

Start with the wn and let say a; > 0. This is bounded by 

\wn{x,k,t)\ < ' ^"W |e-^^-e-^|(l+2/)-™dy. (50) 



/47ri 

Divide the integral to y < \/t and the complement. For the former, Ii , use 

|e 1* — e 4* I < e -»* |e 2t — e 2t I < e |?/| 







to obtain 

h < C\\f\\it-^e-T^e-*''^ < C||/||ii-ie-'=('=^+'^)e-^*'=' 

<C\\f\\it-hMktx)e-i"'" (51) 

where ^ + tk^ > 2kx and k^^ < \/t were used. The integral over y > -\/t is bounded 

,,^L!^r(,e-'^+e-'^)ll + y)-iy. (52, 

V47rt Jvt 

Consider the first term on the RHS. Integrate first over the domain Ix — y | > ■^x. Since 
(1 + y)~™' < t~^{l + ?/)~™^^ and m > 2 we may bound this term by 



-e-^ (53) 



t 

which can be absorbed to (I51> . The integral over |a; — ?/| < i.T in turn is bounded by 

C||/i|ie-''^-'(l + .T)-". (54) 
Since a; > in this domain we may bound i54l as 

CUhe-^'^^ih + hx)-\l + < C||/||ie-''='fc2(l + k^x)-^ 

< C\\fht-ikte~i*''\l + ktx)-"'. 

(55) 

The second term on the RHS of i52\ is bounded by the first one. Thus, altogether, we 
got 

kill < Ct-i {uj{x) + ktuj{ktx))e-i'''\ (56) 



For wi2 we use (I28> . This estimate, for k < t 2 , is readily seen to produce a 
bound like i56i and, for kg > k > t^^ similarly, if we use 
(and replace in i56\ the j by i). 

Thus, altogether, in the domain k < ko, we have proved the decomposition J45t 
with 

|wi(x,k,i)| < Ct-i (lu{x) + ktuj{ktx))e-i"'\ (57) 
For t > 2 and k > ko ( I34> gives 

|(e(*-i)^/)(a;,k)| < C||/||ie-^^-'(*-i)o;(x) < ie-^*e-^(*-i)'='c.(a:) (58) 
upon usmg / e-l^-2'lw(y)dy < Cuj{x). Since 

it can be absorbed into (I47> . 
For i < 2 use 03 and (|70|i: 

Again this and vq can be absorbed into vi . □ 
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Define next w ^ v ~ gi*--^)^ and wi ~ w + vi, i.e., 

V = Vq + Wl ~ Vq + Vl + W. 

We show that w and wi remain bounded in the Banach space B of continuous functions 
s w{s) E Xs with norm 

||w|| :=supsT^||w(-,-,s)|U. (59) 

S>1 

Proposition^follows since v ^ vq + wi. 

Assume now that t > 2. Since the estimates for our semigroup are quite different 
for short and long times, we decompose the integral equation as 

.t-i 

w{t)= Sit-s)auo{vo{s)+wiis)y^ + h,{s)*^)ds + C{t) (60) 



/ S{t~s){^uov{s)*^ + ^v{s)*'^)ds- (61) 
Jt-i 

for later convenience we wrote in the first term v = vq + wi. (We omit the details of 

rt — l 

the case 1 < t < 2. This is simpler: one omits the term with on the right hand 
side of i60\ . and the integral in ( I61l l is replaced by J^. One applies the estimates given 
for C{t) below.) 

We assume that ( I40> holds and that w is in the ball with radius 6 at the origin of B 
and show the RHS of i6Q\ stays in this ball and contracts there. Then ( I39l l follows from 
the contraction mapping principle. Notice that v satisfies by ( I46> . ( I47> and i59\ 

SUp\\vi;;s)\\s<C6 (62) 
s<t 

since |^| < CS. 

Next, recall that for t — s > 1 we have the decomposition 5* = 5*0 + S*! where the 
operator 5*0 annihilates odd functions since Z in ( I25t is even. Since uq is odd and vq 
even, the term involving SoUoVq"^ vanishes and we may rewrite i6Qi as 

w{t) ^ / So{t - s){^uo{vq + wi) *wi + ^v*^)ds 
t~i 

Si{t ~ s)(f uow*2 + ds + C{t) 
= :A{t)+B{t)+C{t). (63) 
We need a lemma about how our norm behaves under convolutions; 
Lemma 6. Let /, g e Xs. Then 

<Cs-^|l/|U||5||. (64) 
Proof. We get from the definition of the norm (using the notation of ( I37> ) 



|(/*ff)(x,k)| < ||/j|,||.9j|, / (c.(x) + |k~p|,c.(|k-pUx)) 

• {uj[x) +PsUj{psX)) 

•(l + |k-p|3s)-"(l+p3s)-"dp. 



(65) 
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Expanding the product the integral gives rise to four contributions. Three of these have 
at least one factor lo{x) and can be estimated by Cui{x)I, where 



/:= y"(l + |k-p|3s)-"(l+p3s)-"dp. 



(66) 



Dividing the integration domain to := {|p — k| < i/c} and the complement Ec, we 
estimate 



I <{l + {\kfsy(^j (l + |k-p|3s)-"dp + j (l+p3s)""dp) 
< C(l + fc^s)-" / (l+p3s)-«dp < C(l + fc^s)-"s-'^. (67) 
The remaining term /' is more compUcated: 

/' := j p,s|k-p|,cc>(p,xV(|k-p|,,a;)(l + |k-p|3s)-"(l+p^s)-"dp 
< Cfc,tj(fc,,x)(l + fc^s)-" / a>(|k-p|,2;)(l + |k-p|3.s)'"dp 

J E 

+ Cu{ksx){l + k^sy j Ps\k-v>\sio{psx){l+p^syA^. 

The first term is of the appropriate form, since the integral is bounded by Cs"'''"^)/^. 
We need to extract a ks also from the second integral. Use |k — pj^ < fcs + Ps- The 
term containing ks is clear, and it remains to estimate 

/ p2^(p,a;)(l+p3s)-"dp< / L+^)'(i+/s)-"dp 

C Gks 



g(d+l)/3 - g(d-l)/3- 

□ 

Let us start bounding the terms in i63\ and consider first C{t). By Lemma|6land 
(I62> we have 

||«*'||.s , ||«*'||.s < CS^s-^. (68) 

Thus the bound i36\ gives 

\C{t)\<CS^j' dse-^''"^'-'\t- s)~i y"dye-(*-^)"*l"-^l 

• (uiy) + ksLuiksy)){l + sfc3)-"s-^ . (69) 
The following simple estimate is needed repeatedly below: 
Lemma 7. Letk,p> 0. Then 

f , _i.|t-7/I / N ^ \ Ck'^uj(px) fork>p, ^^^^ 

/ dye ''I'' ^^Lo{py) < < ,_ , , (70) 
J I C{p uj{kx) + k u!{px)) for k < p. 



12 



Proof. For k > pwe decompose the integral to the set E given by |a; — y| < ^\x\ and 
its complement Ec'. 

J Je J Ea 

and ^ follows. □ 

Using Lemma0and e"'*"*-' ' 1^"^' < e"'^"''' we may now bound the y integral 
in ( I69l l by 

dye-l^-«l(w(2/) + ksuj{ksy)) < C{uj{x) + ksUj{ksx)) 



Since s E [t ~ l,t],hy changing the constant C we may replace s by t except in the 
t — s factors and end up with 

\C{t)\ < C5h-^ {uj{x) + ktuo{ktx)){l + tP)"" 
i.e. the Bt norm is bounded by 

\\C{t)\\t < C5h-^ (71) 



For A and B in ( I63> we need to distinguish between the various k values, 
(a). Let first k <t^^. Using (|68} and (|3T} we bound B{t) by 



\B{t)\<CS'{h+L2+l3) 



t-i 



h / ds / dy {t - .s)-'e-^\^-yls'^ (uiy) + k^k^y)) (72) 



d-l 



t-1 



h ■■= J dsj dy {t - s)-'e-'VT^-'^y^s-^{uj{y) + kMhy)) (73) 

h ^ dsj dy {t - sy^^^e-'^s-^{u;{y) + k^k^y)) (74) 

For Ii we use Lemma0 s~'^l'^ < kg < 2.s~^/^ and uj{x/y/s) < uj{x/ \/t) to deduce 

h<C I ds{t~ s)-^s-^{u{x) + ^Lo{^)] 
Ji ^ Vs / 

<£}^L{x) + kMktx)). (75) 
i 3 V / 

For I2 we use simply 

Ji {t-s)s— Vt t— Vt 

assuming d < A. (The bound is not useful in the case d ~ 2.) For all d > 3 we obtain 

h < Cr^/^\ogtktUj{ktx). 
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For I3 we use Lemma0again, then estimate both Lu{x/y/s) and uj{x/\/t — s) by 
uj{x/\/t): 



<^(-(x)logt + .(^))+^.(^)) (76) 

for d < 4. Again for d — 2 the resuh is not useful. For d > 4 one obtains I3 < 

Ct-^{uj{x) +uj{x/Vi))- Thus for d > 3, 

\B{t)\ < CSh-'^"^ log* (^Lo{x) + t-^Lu{^)^ . (77) 
For the A{t) in (I63> we use Lemma|6lto bound 

and then ( I30> to get the bound 

• s""^"* (w(y) + k,Lu{ksy)) 



which for d > 3 is bounded by 

CSh-^uj{x). (78) 

Thus for k < t^i 



\A{t) + B{t)\ < CS^ \ogt + . 



(79) 



(b)Let us next consider k E [t 2 ^ ko] and start again with the B{t) in ( I63> . Using (I33> 
is bounded by the sum of 



C6^ I ds / dy ((i-s)-ie-^l^-«l + (t-s)-^e-"'^l^l-"l^'l + (t-s)-='/2g-cfe|^-j/|^ 

^ + fc,w(fc,y))e-4(*-")'='(l + s/fc^)"" (80) 



S 



and 

• (a;(y) + fc,t^(fc,y))e-3(*-^)^' (1 + sfc^)"". (81) 

The first integral resembles the case k < t^^/^. To estimate (I80> use 

e-5(*-')'^''(l + sfc3)-" < C{l + tk^)-", 
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and assume d < 4. Then (18 0> is bounded by 

t-i 

C6^n + tk^y [ ds— (u}(x)+uj(kx)+ / (uj(x) + ksUj(ksx))) 

J s—{t-s)^ k^t-s J 

< CS\l + te)-^(^^{u;{x)+Lj{ktx)) + -^{VtLo{x)+Lu{ktx))). (82) 
In the case d > 4 replace (d — l)/3 by 1. For d > 3 we thus get the bound 

Cb''(\ + ifc3)-»^(^(a;) + kt^{ktx)). (83) 
Thus (18 U remains. Using ( I70t we have 

y e-5fc|^-3/lw(2/)d?/ < C7(A;-iw(a;) +w(fca;)) 

and 

y" e-^'^l^-''lfc,a'(fc,?/)d?; < C7(fc"ifc,tj(fc,a;) + ^(fca;)) 

whereby 



Ay (fc^e-'^'^l^-^'l + fc3e-=('=l^l+l^l) + fc^g-sl^^-!^!) (^(y) + ksUj{k,y)) 

< Ck^{uj{x) + ksu{ksx) + kuj{kx)) < Ck^{uj{x) + {kt + s"^)Lj(fcta;)) 
since Lo(kx) and Lo{ksx) are bounded by Cu){ktx) for /c > t^^/-^. Thus 

/t-i 
d.s s-^e-3(*-^)'=' (1 + sk^ye {lu{x) + (fct + s-^)u;(fcta;)) . 

Using 

e-i(t--)fe'(l + sfc3)-"A:2 < C{t - + tfc^)-" 

we get 

/i-l . 
dss-'^(t-s)-2/3(^l + (-)^j. 

The integral is, for d = 3, bounded by Ci^ s and we end up with the bound for k S 

(ED < C(52t"^(w(a;) + fctw(fcta;))(l (84) 
A{t) in (I63> is for fee [^^^ , ^o] bounded by the sum of 



+ fc,w(fc32;))e-3(*-^)'^-'(l + sfc^)-" (85) 



s 
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and 



CS^ J dyJl 'd.s((t-5)-5e-5(l-l+l?'l) + (t-.5)-ie-3l-l--'=|j'l) 

In the former expression we replace as above fc by (t — s)^^^^ and end up with the 
bound 

Cdh-'^/^^uj{x){l + tk^)-''\ (86) 
and the latter is smaller with <~^/^^ replaced by t"^/^. Thus for fc € [t~K ko] 

\A{t) + B{t)\ < CSh-^/^^iujix) + ktuj{ktx))il + tk^y. (87) 

(c) Finally we suppose k > ko and use the bound ( I34> and Lj{ksx) < Cuj{x) which 
holds for such k to get 

/t-i j- 
ds J dye~^^'(*-^)e-=l"-^ls-'^w(2/)(l + sfc3)-" 

< CSh-^/^uj{x)il + tk^y (88) 

The bounds (|7D, and (|88j show that our ball is mapped to itself. The con- 

tractive property goes similarly. 

4 The spectrum 

In the remaining sections we will present the linear analysis needed for the semigroup 
kernel estimates in Lemmasm2land|3lof Section|2l Since the linearized Cahn-HilUard 
operator is of the fourth order, the calculations are lengthy, and we will be brief in the 
most tedious details. The full details may be found in the thesis 0|. The linear analysis 
is based on an approach outlined in |6J. 

We remark that, in addition to the existing notations, boldface letters will mainly 
refer to vectors in and 4x4 matrices in the sequel. A boldface symbol and the cor- 
responding italic symbol are not always mathematically related: consider for example 
R and R, or v and v. 

The following lemma summarizes the properties of the spectrum of DkH^ alluded 
to in Section|2] 

Lemma 8. For any k > Q the spectrum of DkHk is contained in [fc^, oo). There exists 
an e > such that for < k < e the bottom of the spectrum is an isolated eigenvalue 
Co of multiplicity one with Co < Ck^. The remaining part of the spectrum is bounded 
from below by jk^. 

Proof. DkHk has asymptotically constant coefficients, thus the essential spectrum 
is [fc^ + fc'*,oo), the same as that of the constant coefficient limit + Dk, see 
lis] Proposition 26.2]. Eigenvalues on the other hand are the same as those of the 
self-adjoint operator Ak := d]/"^ HkD]/"^: the eigenvectors are mapped onto each 
other by Df!", which is invertible when fc > 0: we have the convolution kernel 
{Dl^){x) =:Gfc(x) = (2fc)-ie-^-|-l. 
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Hq has two isolated eigenvalues at and 3/4 and a continuous spectrum [1, oo), see 
||9] page 79]. The eigenfunctions are V from il6i and x i-^ sinh(a;/2) / cosh(a;/2)^. 
The spectrum of Ak can be studied with the minimax principle: 

. „ {u, Hku) 3 
Co = mf , , < , , < Ck^ 

u {u,Gku) {V,GkV) 

for small k. We also get the lower bound: 

« (it, fc^Gfeu) 

To see that Co is a discrete eigenvalue of multiplicity one we proceed further with the 
minimax principle: 

. {u,Aku) . {u,Hku) 

Ci = sup mi — — = sup mt - — — — - 

V uev^ {u,u} V {u,Gku} 

(Gf,' u,v)=0 



ju inf (^^^ku) ^ ^ ^ .^^^ {u,Hou) k^\\uf ^ 

V tiGti-L {u,GkU) ~ V iiGu-L {u,k'^Gku) ~ ^ 



which for small k is larger than Cq. Here we used knowledge of the spectrum of i/o: it 
has zero and 3/4 as isolated eigenvalues and a continuous spectrum [1, 00). □ 



5 The resolvent: a first order system 



In Sections ISHU we analyze the integral kernel of the resolvent of the operator C, — 
DkHk, where C is a complex number outside the spectrum of DkHk', see (I20> . We 
start by writing {C, — DkHk)u = / as a system of ordinary linear differential equations 



u'(a;) = A(C, k, a;)u(.T) + b(a;), 

where 

M^{u,u',u",u"'f, 1 

and (for V, see Mbl ) 



(0,0,0,-/)^ 



(89) 



(90) 



/ 



Vc 



10 0^ 

10 

1 

^-k'^V + V" 2V' l + 2k^ + V Oj 



(91) 



We define the x-independent matrices Aoo(C,A:) :— lim2;__>±oo A(C, fc, x), and 
R(C, fc, x) := A{C, k, x) — Aoo(Cj the latter can be bounded by Ce"'^' for k) 
in any compact set. The positive constant G may depend on the set. 

The eigenvalues of Aqo are the zeros of the polynomial 



(:^{-^,' + k'Y-{-^,' + k') 



which are 



(92) 
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Unfortunately they have somewhat poor analyticity for k and ^ near zero. This problem 
can be overcome by writing 



(^k^ + k* + {l + 2k^fT^ 
and using A = (fc, r) as the parameters. Then (19 1> and \92\ become 



10 0^ 

10 

1 

+ 2fc2)2r2 - fc2^ + V" 2V' 1 + 2fc2 + 1/ 0/ 



(93) 



and 



To fix the branches when the second sign is negative, set 



(94) 



± 



4^2 



±ir 



1 + 1^1+47^ 



(95) 



and choose 112 to have the positive sign with the principal branch in the denominator. 
Now the are analytic functions of k and r when fc e C \ ±i[l/V2, 00) and r G C \ 
±i[|, 00). We are mostly interested in real k. Then the cuts in the t plane correspond 
to C € (—00, —1/4], values which we do not need for integrating around the spectrum 
of DkHk- T e ±i[0, i/c] U M is mapped onto [3fc2/4 - k*^, 00) which by LemmaEl 
contains all of the spectrum of DkHk except for the lowest eigenvalue, which should 
be somewhere around r w i/c when k is small. For larger k it is more useful to use 
the A:^ bound: t G ±i[0, k/ (1 + 2k'^)] U M is mapped onto [fc*, 00) which contains the 
entire spectrum. 

The number r is mapped two-to-one on ^ and we will mostly keep r in the upper 
half plane. There Re /i2 < when k is real. However, we need some results to be valid 
also in a small complex neighborhood of the origin for both k and r, because we will 
develop a power series there later 

The number fii will be the eigenvalue with the most negative real part, /i4 := — /ii 
and ^3 := — /Li2. We have 



M?4 



Ijl^l + 2k^, fiifj,2 = -iT{l + 2k^), 
- i^lI = (1 + 2A:2)^i + 4r2 = ^1 + AC,. 



(96) 



We denote by Vj and Wj the right and left eigenvectors of A^o'- AooVj = l^-j^j 
and Wj Aoo = l^j'^j- The vectors Vj and Wj are easy to express in terms of /ij and are 
thus also analytic. Normalize the eigenvectors so that (vj)i = (wj)4 = 1. This results 
in 



Vj = (l, Mj, Mj, f^i 
w, = {11,(^1^-1-2^), 



1 - 2k'\ 
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6 The resolvent: solutions of the homogeneous equa- 
tions 



The resolvent for (I89> will be constructed using the solutions of the corresponding 
homogeneous equation and its transposed equation (A(x) :~ A{(, k, x) and so on): 



So, this section contains a study of the solutions of ( I97> and ( I98l l. The motivation for 
(I98> is that if u and z are solutions of ( I97> and ( I98> the product z(a;)u(a;) is independent 
of X. 

As the coefficient matrix A(.t) tends to the constant Aoo when a; — ^ oo, one would 
also expect the solutions of ( I97t and ( I98> to tend to the solutions of the corresponding 
constant coefficient equations. Indeed, Pego and Weinstein show in f6i\ that if there is 
a simple eigenvalue with a smaller real part than the other eigenvalues, ( I97t admits 
a unique solution with the property that e^^'^^u^ (x) ^ vi as a; ^ oo. Similarly 
( I98> has an unique solution with e^i^Zj^(x) ^ wi as x ^ — oo. 

We now examine solutions corresponding to the remaining eigenvalues yUj for j > 
2. We prove things for but the z~ case is similar. Fix a set of values A :~ {k, r) C 

to work with. (As explained in the previous section, the dependence of A or Hj on 
A or A: or r is analytic in the natural way, if A is restricted to a suitable domain.) Choose 
some A-dependent n and substitute u(a;) =: e^^\'{x) in (|^} to get 



where B := Aoo — M ™d Il{x) := A(a;) — Aoo- 

We want to divide the eigenvalues of Aoo into two sets: those with smaller real part 
than /i and the rest. However, we need some margins in order to keep things uniform in 
A, because the eigenvalues may not maintain the same order for all A (see Sections- 
Let Ei := Ui^i '^*5r(/ii — Aoo)". Choose constants a and (3 and a set / C 
{1, 2, 3, 4} so that Re(/ii — ^) < a < when i e / and Re(/^i — IJ) > P > —p 
whenever % ^ 1. The numbers a and fi are assumed independent of A. Let P be the 
projection onto ® Ei which commutes with Aoo and let Q := 1 — P. Fix some 
positive xq and define 



for bounded continuous functions v : [xq, oo] ^ C^. The proof of the following facts 
is straightforward: 

Lemma 9. For sufficiently large xq, T is a contraction in the norm 



v^!{x) — A(a;)u(a;) 
z'(x) = — z(.T)A(a;). 



(97) 
(98) 



v'(x) = (B + R(x))v(x) 



(99) 




(100) 



sup |v(x)|. 



Corollary 1. {Tw){x) = 0{e"'^'''^-P'°'+'^'=) for large x and any e > 0. 
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Hence, v = v + JFv can be solved for v given any v. For such a solution 

(v - v)' = PRv + BTpv + QRv + BTqv = B(v - v) + Rw. 

In particular if v is a bounded solution of v' = Bv (a constant coefficient equation) 
for large x then v will be a solution of i99\ with the same asymptotic behavior as v in 
the sense that v — v = tends to zero exponentially fast as a; ^ oo. 

Corollary 2. For each eigenvalue fii there is a solution of \91\ which behaves asymp- 
totically (as X — > ooj like e^"^v.i. 

Proof. Pick /I = yUi, a = — |, /? = — |, e = and v(a;) = v^. It is clear that a 
suitable / can be found and can be used to get a solution for x > x^, which can then 
be extended to the whole real line. □ 

In general the solutions of Corollary |2l are not unique even after fixing normaliza- 
tion of V, because one can add similar solutions corresponding to /i^ with smaller real 
parts than /i^. 

Corollary 3. If the and Vi are analytic functions of \ in some domain ofC^ and we 
can fix the I in the preceding proof uniformly for all X, then the solutions in Corollary^ 
are analytic in this domain when evaluated at some fixed x. If Mi is such a solution then 
|e~'''^'^^^Ui(A; a;)— Vi(A)| < Ce^^^ for x > and X in compact subsets of the domain 
( C depends on the subset). 

By Corollary 121 for j £ {1,2,4} there is a which solves the homogeneous 
equation dxUj' = Au^ and behaves like e^'^^v^ as x — > oo. There is also a zj 
which solves the transposed equation dxzj = —zJ A and behaves like e^^^^Wj as 
X — oo. It would also be possible to define U;^ and z^ in this way but that would not 
be very useful because Aoo is defective at r = 0: the eigenvectors V2 and V3 collide 
and we would not have a set of four linearly independent solutions. Thus with some 
abuse of notation we require instead 

uT(x) ^ and Zn (x) ^ 

^ - /i2 1^3-1^2 

as a; — ^ cx) and x —00, respectively. As ji^ — /i2 ^ these converge (pointwise in 
x) to solutions with linear asymptotes. 

Our expression for the integral kernel of — DkHk)^^ will contain only and 
z J for j E {1,2}. However, the other two values of j are needed for understanding the 
behavior of the kernel. 

We summarize the properties of the solutions of the homogeneous equations in the 
following theorem. 

Theorem 2. 

Rc(^i±^2)<-| and Rc(/^2) < ^- (101) 

Then Wt\ with ( l9ll has solutions uf such that |u^(a;) — e'^'^Vij < C\e'^'^'~^^^\ for 
i G {1,2,4} and 

|u^(x) - ^ 1 < C|e(^-^).| 
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when X is bounded from below. Furthermore, 

|9>+(x)-/.re^'"v,| <C|e(^'-^)-|, ne{l,2}, 
|5,u+(x) - i(e^^^V3 + e"^^V2)| < C|e(^-^-^)-|, 

Similarly, (|98j /lai solutions z~ iMc/i f/iaf |z~(a::) — e^^*^Wi| < Ce^^~^'^^ for i e 
{1,2,4} anc/ 

w/ien X is bounded from above. The corresponding estimates for the derivatives hold. 

For i S {1,2,4} the functions u,^(x) and (x) are analytic in A wherever the 
assumptions hold. The functions {x) and (x) are continuous and when 

Re(//2) > (102) 

also analytic. The constant C above depends on X but can be fixed in any compact 
subset. 

The proof is basically just application of Corollaries|2land|3]with some additional 
details for the z = 3 case. There we need to use different choices for the set / depending 
whether Re(^2) is near zero or not. The results can be glued together with a partition 
of the unity but analyticity is then lost. See for more details. 

Specifying the assumptions in terms of /ii and /i2 was perhaps a bit opaque. How- 
ever, there are two cases which interest us. The first is A « 0. In a small enough 
neighborhood of the origin both dlOlt and J102> clearly hold. The second case is when 
k is real and r in the upper half plane. Then /i2 is in the left half plane. For the first 
inequality of dlOU use ( I96> to get (/ii ± 112Y = (1 + 2A:^)(1 =F 2iT). From this we see 
that the inequality is equivalent to r lying beneath an upwards opening parabola with 
apogee at i(128fc-^ + 39)/(256A:^ + 128). This is always above the values corresponding 
to the spectrum. 

To end this section we still note the following fact, which is a consequence of the 
symmetry of the Cahn-Hilliard equation under the reflection x ^ —x: 

Lemma 10. Under the assumptions JlOU there are also solutions uj and z'j' of i9H 

and ( I98t . respectively, with prescribed behavior for x —00 and x 00: (e.g. 
u^(a::) ^ e^^^^V4 as x —00). These can be expanded as linear combinations of 
and zj. The coefficients of this expansion are continuous and, when ( I102t holds, 
also analytic. 



7 The resolvent: a resolvent formula 

We now want to write solutions of ( I89l l. that is, of u' — Au + b, using the solutions 
of the homogeneous equations. We will use the following resolvent formula: 

u{x)= r u+{x)in+)-'z-iy)h{y)dy+ P {x){n-)-'z+ {y)h{y) dy, 

J —00 J 00 

(103) 

where = (uj*^, uJ), = (^0 ~ z^u*. These products are indepen- 

dent of X and it is also easy to see that fi^ = — For simplicity we will from now 
on denote fi+ just by ^l. For ( I103t to make sense il of course needs to be invertible. 
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Theorem 3. Asswne Re fii < Re fj.2 < ond C G p(DkHk) ( the resolvent set). Then 
= z~u+ is invertible. 

Proof. Assume the contrary: let Vl{ai,a2)'^ = 0, i.e., z~u = for i E {1,2}, 
u := aiu^ + a2uj ^ 0. The function u can also be written as J2'j=i l^j'^J ■ have 

I Vi otherwise 

as X — oo. Thus {x)uj (x) (54jWiVi but the left side is actually independent 
of a;. Since /ii is a simple eigenvalue we must have wivi ^ 0, consequently Zj^u = 
implies P4 = 0. In a similar vein z^u = implies (3^ — 0. Hence u = /3iUj^ + 132^2 
decreases exponentially at ±00, giving a nontrivial solution to {( — DkHk)u — 0. 
Such a solution was assumed not to exist. □ 

We omit the proof that ( I103> is actually a solution of ( I89l l. 

The solutions can be written in terms of a scalar valued function . A straight- 
forward computation yields: 

Lemma 11. 

u± = (C/±, d^Uf, dluf, dluff 

zf ^{{Hkd, + Pd,-v')zj, -{Hk + e)zf, -d^zf, zf) 



with (C - DkHk)Uf = and (C - HkDk)Zj = 0. For j e {1,2,4} we have 
lim,^ooe-''^"[/+(.T) = lim,__ooe^^^Z-(a;) = 1 while \im,_,oo e^f'^'^U+ix) ~ 

^12"^ ~ liiiix^-oo e''^^Z^(x) — - — !^^^ = {when 112 — 0, replace the fractions 
by their limits, i.e., x and —x). 

Recalling ( I90l l we get for the original equation: 

Theorem 4. Under the assumptions of Theorem^the integral kernel of the resolvent 
is given by 

Rix,y) := {C-D,H,r\x,y) = l-U^(-)^-'Z-iy) Mv < ^, 
^ ^ i^J y ^yj \-U+{-x)n-^Z-{-y) fory>x, 

where U+ = {U^ , U^) and Z- = {Z^.Z^^V- 

Note that from {x)Vl^^z^ {x) — vT {x)Vl^^z'^ {x) = 1 it follows in particular 
that the (i, j) component of the left hand side is for i < j. From this we see that our 
resolvent kernel has continuous derivatives with respect to 2; or y up to total order two. 



8 Estimates for the resolvent 

The leading terms of the semigroup will arise from the resolvent with small k and |^| 
values. Deriving estimates for the resolvent kernel ( I104> for these parameter values 
is the main task of this section, leading to Lemma [T2I and Theorem |5] After this we 
consider large parameter values briefly in Theorem|6l 



22 



We use ( I93> and assume in the following that fc < e, |t| < e for a small enough e. 
In particular we assume that the smallest eigenvalue is isolated (see Lemma|8j and that 
the analyticity condition (I102> and the other assumptions of Theorem |2] are satisfied. 
We will proceed to develop 57 of (I104> into a power series to get some explicit leading 
terms for the resolvent. Thus this section is mostly about computing derivatives of 
things at fc = r = 0. 

We denote the solutions of the homogeneous equation at fc = r = with a ring. 

p(i = — 1, /i2 = and 



U+{x) 



1 



4cosli(f)2 



o , -l-6e^ + 5e2^ + 2e3^ + 6e2^a; 
U2 (^) = 9..n+.-^2 (105) 



2e^(l + e^) 
Zr(a;)=log(e- + l) 
Z^{x) = 1. 

By a mechanical computation we get 



1 




The solutions which are asymptotically equal to x i-^- |x|, i.e., U^{x)/x 1 as 
X oo and {x)/x — 1 as a; ^ — oo can also be computed: Z^{x) = ~x but 
the expression for is lengthy and we omit it. The following products will be needed 
later: 



■■'1 



z 



2 




u+ = ( -1 I and Z3" (u+ u+ u+) = (O -1 O) . (106) 
We will also also need the rapidly growing solution lJ^{x) — 4cosh(|)2. 



First order: Define v^(j;) = e~^^^\i^{x) and wj {x) = e^^^^zj (x) for j e {1, 2}. 
By Corollary |3] they are analytic functions of A for each fixed x and satisfy |v^(a;) — 
Vj I < Ce^^/2 for X > and | w~ (x) — Wj | < Ce^/^ for x < where C can be fixed 
independently of fc and r. They will also satisfy the obvious differential equations, 
which we now differentiate, then set fc = t = 0: 

dxd\v^ = (A - fij)dxv^ - dxfj.jv^ 

(note that 9a A = 0). These are just inhomogeneous versions of the equations satisfied 
by and wj . The "initial conditions" are again asymptotic: \dx'v^{x) — dx'Vj\ < 
(j^-~x/2 ^jjj \dx'wj (x) ~ dx'Wj I < Ce^/2 by Cauchy's estimates. The first component 
of dxVj and the last component of 9a Wj are zero because of the normalization used. 

Using all available information developed until now, we end up with these reason- 
ably simple expressions: 

dxzj = dxfijz^ (107) 



9au+ = 9am,u^ 
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and, using (I106> . 

dxfl = 9az u+ + d\u^ = 
with dxi-ii = dkii2 = and 9r/i2 = i- 





-2^x^Ji2)' 



Second order: Above we obtained non-zero leading terms for the right column of Vl. 
We go on with the left column in order to eventually compute some non-zero terms for 
dot Vl. We have to deal with 

since 9a = at A = 0. It is convenient to estimate the terms on the right hand side 
at separate values of x, which is possible with the following trick: 

dM^-xxt) = -dl{z-A)n+ + dlz-Ant = -z^dlAuf 

(recalling d\A = 0). Now we can write 

dlizrnt) = idlzrnt)ixo) + {zV dlnt){x^) - T z-9^Au+ dx. 

Jxa 

Taking the limits xq ^ —oo and xi ^ oo and using the asymptotic behavior of Uj'" 
and other functions involved everything except the integral vanishes and we get 

dl{z-n+) = , dA{z-n+) = 0, dl{z-n+) = 

Third order: Similarly we obtain S^Zj^Uj*") = 0. 
Collecting everything from above we get 

Lemma 12. 

det r2 = fc2 + + 2iT^ + |irfc2 + 0{k'^) + 0{t^). (108) 

According to (I104> . R{x, y) = R{-~x, — y). Thus it suffices to consider the case 
y < x. Then 

R{x,y) = ^U+ix)l5Z~{y). 

where 

/ -2ir -1 \ /0(A2) 0(A2) 

^-[P+T^ -IP - ry ^ [oiX^) 0(A3) 

In our small neighborhood (fc, t) « the only possible singularity is that det 51 may 
become zero, producing a pole in the resolvent. To examine the x and y dependence 
of the resolvent it is easier to work with F{X; x, y) :— —U^{x)l3Z^ (y) which has no 
such singularity. 

FiX;x,y) = U+ (A; x) {Z^ (A; y) + 2iTZ^ (A; y)) 

2 2 

i=\ j=i 
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where the 7^ are analytic functions of A = (fc, t) and 0{\^). It turns out that we need 
to look at the U2 term a bit more carefully so let us write it more explicitly: 

F = U^Z^ + TitU^Z^ - (fc2 + + 0{X'^))U^Z^ 

2 

+ E ^iJ^i^i + ^22U+Z^ . (109) 

As these expressions for F get longer we drop the arguments to reduce clutter U is 
always evaluated at x, Z aiy and everything depends on A. Also, U~ {x) := U^{—x) 
and Z+{y) := Z-{-y). 

When y < < a; it is easy to bound ( I109> by using Theorem|2]but in the other cases 
where y < x we do not know much about the behavior of either (x) or Zj' {y). We 
can get around this by using Lemma [Tol 

When < y < X write 

4 

Z-{X;y) = J2h,{X)Z+{X-y). 

3 = 1 

Explicit computation yields 

hij - 5ij + + 0{X) 

b2j = S2j - 2ir<53j +0{X^). 

Thus 

2 4 

F = U+Z+ + 2iTU+Zt + E E I^^^UtZt 

1=1 j=i 

for some coefficients /Sij, which are analytic functions of A and 0{X^). F must be a 
bounded function at least in the resolvent set. Hence /324{X) must vanish for all A. The 
j = 3 terms, however, are a bit tricky: we have 

|e-'^3.^3+(y) - — 1 < Ce-^y 

and /i2 ~ ^T. (323 = — + 0{X^), which is more important than /3i3 because 
of the rapid decay of . Let us write F a bit more explicitly and also bring in the 
U^Z+ term: 

F = U+Z+ + 2iT[/+Z+ - (fc2 + + 0(A=^))(C/2+^i+ + U:^Z+) 

4 

+ E /^ij^i^^.^ + P22U+Z+. (1 10) 

When y < a; < we do the same thing but with U: 

4 

aij =(5i, +0(A2), 

a2j = 8(5ij - 52j ~ \5i.j + 0(A). 
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Plugging this in we get 

4 2 

with some coefficients (A) = 0{X^). Again Q!42 must vanish to keep the kernel 
bounded. To compute the Z2 term coefficient to second order it seems that we 
would need 013 to second order, which seems difficult to compute. However, we can 
get around this difficulty by using twice continuous differentiability of F when y w 
X <C 0. The limits 

lim F{x,y)=a:i2U^{x)Z-{x)+aiiU^{x)Z-{x)+0{e't''\''\) 

lim F{x,y) ~ lim F{—x,—y) 

= P23U+{-x)Z+{-x) + /3i4U+{-x)Zt{-x) + ©(e^^'^l^l). 

must be equal. Working in the region where Re /j,2 < and taking the Hmit a; ^ — 00 
yields 

"32 , /323 , a 

h a4i = h /?4i- 

2/X3 2/X3 

Repeat the same for d^F: 
lim dlF{x,y) = ai2dlU^{x)Z:,(x)^a^xdlUl{x)Z^(x)^0{e^^^\^\) 
lim d^F^x, y) = lim dlF{~x, —y) 

y — >x+ y — >x+ 

= (32sdlU+{-x)Z+{^x) + pudlUt{-x)Zt{-x) + 0(e2A'2|.|)^ 



hence 



The two equations for 032 and a^i are linearly independent (recall \9()V ) and their 
obvious solution is aij ~ /3ji. Thus we get 

F = U^Z^ + 2iTL/f - (fc2 + + OiX^)){U;^Z^ - U^Z^) 

2 

4e{l,3,4} i=l 

To combine the expressions for F in different regions we denote by fij{x, y) the 
part containing Uf^ i^)^f iv) < I109> . ( II lOt or ( II 1 1> . We shall also get rid of and 
Z~ by using absolute values. We write 

4 4 

with, e.g., /n(x-,y) = {2iT + 0{\^))U+ {\x\)Z+ i\y\). 

Thus we have written F in terms of and Z^^. For small A these functions can 
be approximated by their explicit forms at A = 0, given by ( I105> : 
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Lemma 13. 



Proof. By Theorem|2| 



Z+{y) = e(^i+i)^Z+(y) + 0{\^e^^''-^^y) 



for i S {1, 2} and similarly for the . Recall that Z2{x) = 1. In some cases the 
first A-derivatives of e^'''^[/,^(x) and e^^^'^ Z^ {y) also vanish at A = 0, hence the 
remainders are 0{\^) rather than 0{X). 

By TheoremlJl 



□ 



U+{3 



similarly for Z^. We also need a couple of derivatives with respect to y: 
Lemma 1 4. For n e { 1 , 2} 

\d^z+-d^z+\<c\x'e-^y\, 



Proof. e^^'^ydyZ2{y) — known to be ©(e"^/^) by Theorem|2]and vanishes at 
A = 0. Its A-derivative at A = can be computed (we computed d\Z2 earlier) and 
also vanishes. Z^ is similar except without the A-derivative. For Z^ use 



□ 



To simplify further we use /i2 = —^3 = iT(l + ©(A^)). Let c be an upper bound 
for the O(A^) term, assumed to be conveniently small. If r is in the sector {Im(T) > 
2c|r|} and y > we have, denoting z :~ vry, 

|gP2y _ giryj ^ |e-(e'^(^')- „ i)| < CjA^e^+^l^l | < C\\^e^'^y\. 
If in addition x — y> ^ 



2^2 



2ir 





< c 








T 
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Using these results we now estimate dyfij for the of and n E {0, 1, 2}. 
In the case ?; < < a; we have from ( I109> by Lemmas[T3land[T4lthe following terms 
and estimates: 

dyfi2 = ?7i+(a;)a;e'^l^yl +C'(AV"e-^l^l+^'^l^yl) +C'(A2e-5(l^l+|y|)), 
a,"/n = 2iTU+{x)dl^Z+{\y\) + ©(A^e-^d-l+l^'D), 

ay"/2i = -(fc' + T^ +0(A3))(e'^l-la;'i+(|y|) (113) 

+ 0(e-5(l^l+lal)) + ©(A^es'^l^l-sl?/!)) , 

5;'/22 = ©(AVe^'^d^l+lfl)). 
When < y < a; equation (II 10> produces equations (II 13> plus the following terms: 



2ir 

+ C'(AV"-^e5'^l^-yl) + 0(6^51^-2/1) 
+ C'(Ae^'^l^-^l-5l2'l) 

When ?/ < a; < equation dl 1 W gives the terms of dl 13t except for a sign change in 
/21 and 

a^"/3i=0((|fc2+r2| + |A|3)e-il-^l), 

/ piT(|a;| + |y|) piT\x~y\ 

d-hi - -(fc' + + 0{\^)) [d^^ ^— 

+ ©(AV-^e^'^l^-yl) + ©(A^e-^l^-^l) 



9;;/4i =0((|fc2+r2| + |A|3)e-3l--^l). 



All the other /^j are zero. The case x < y reduces to the preceding cases by F{x, y) = 
F{—x, — y). We treat the leading terms separately and summarize: 

Theorem 5. There exists e > such that for |A| < e the kernel F, where F :— 
(det n)R = Fn + Fi and Fi = Fiq + Fn + F12, satisfies 

Fo{x,y) ^ U+{x){e^^\y\ +2irZ+{\y\)). 

' ~^^ie'^^''+y^ - gi^k-yl) when xy > 0, 
when xy < 0, 

F,i{x, y) = -sgn(x(x - y))(fc2 + T^)e'-\-\ Zt {\y\) 
\d'^Fi2{.x,y)\ < C(|A3T"-ie5'"l^-2'l| + jA^e-^l^-^'l | + lA^e^'^l^l-^l?/! |) . 

Fq has two continuous derivatives with respect to y. Although we have tried to 
write these expressions so that they would be valid everywhere the other pieces are 



^"10(2;,^) 
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a bit rough around the edges. Hence the estimate for the derivatives is only valid 
for X ^ y ^ {) ^ X. F itself has two continuous derivatives, see the remark after 
Theorem|3 

We finish this section with an estimate for R := {C, — DkHk)^^ when |C| is large: 

Theorem 6. Let k < e. For any a € (0, ^) there exist C, c, r > such that for all t 
with \t\ > r, argr £ (a, n — a) and n S {0, 1} 

\iRD]:)ix,y)\ < -^^e-^VF\\^-y\. 

\t\2 " 

Proof. We use the second Neumann series 

oo 

(C - DkHkr' ^{C-A- B)-' = ^((C - A)-'By{C - A)-^ (114) 

i=o 

with A -.= 01 + Dk and B := DkV. Let i?oo (C - A)-\ \V{x)\ < 6e-l^l but we 
also need exponential estimates for the kernels of i?oo and R^oDk- It is easier to work 
with their Fourier transforms. The poles of R^o are i/Xj where the are given by \9A\ . 
Expanding in partial fractions we have 

- , . I f I 1 



RooDkip) 



When |r| is large j = 0{^/t) and by \96\ the coefficients of the partial fraction expan- 
sions are O ( 1 /r ) and 0(1) respectively. Assume | Rc /x j | > 2a for some a. Then 



/ 



1 dp t: 



Thus \RaaDk{;x, y)\ < ^e""!^"^!. If r is in an appropriate sector away from the real 
axis we can choose a proportional to -y/r and get convergence for the Neumann series. 
The estimate of the theorem then follows easily. For slightly more details see 0. □ 



9 Estimates for the semigroup 

To complete our paper we are left with the task of deriving the semigroup estimates of 
Section|2j mainly from the resolvent estimates of Section|8l However, when t is small 
or k large a crude estimate, following from standard Fourier analysis and perturbation 
theory, will suffice: 

Theorem 7. There exist C,c > such that for any a £ (0, 1] 

C 



|(e-*^^^^Dfe)(a;,y)| < _^e("('^^'^)+^)*-''l"-yl 

where t > 0, x, y G R, and 

a{k, a) := -^Uk^ + k^) + c{a^ + a^). 
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Notice that for bounded t we can choose a = t ^Z** to get rapid decrease in |a; — y|, 
or for large k we can choose a = ek for an e such that a{k, a) + c < — |fc^. 

Proof. Another perturbation theory argument with A = + Dk and B = DkV. This 
time use the expansion 



n=0 



n=0 

r„(i) = (-1)" j ■ ■ ■ j g-(t-t.)AB^-(t^-t„^i)A^ . . . g-tiA . . . j^^^ 

{0<ti<-<t„<t} 

Fourier transform and the usual imaginary translation trick yield 

t4+ 2 

with some C,c, M > valid for all a > 0, j G {0,1}. By induction we get 

C(2Afr(i))" 



V 4 ^ 



The prefactors can be estimated by the power series of e'^*, yielding the theorem. For 
more details, see |4l. □ 



For other parameter values we use the Dunford-Cauchy integral i20\ and the resol- 
vent estimates. There we have more cases: for small k and C, we use Theorem|5l for 
large ( Theorem|6land for other k and C just Theorem^] 

Let e > be as in Section |8] and Theorem |5] assume < A: < e and that t is 
bounded away from 0. There is a spectral gap between the lowest eigenvalue at 0{k^) 
and the rest of the spectrum from 0{k'^) upwards. By analyticity, the integration path 
in ( I20> can be modified to consist of two parts: one surrounding the smallest eigenvalue 
and yielding a term which we denote by A'P"'*^; another part surrounding the rest of the 
spectrum and yielding a term X"''". 

We first estimate K'^°^'^. Changing the integration variable to the r of ( I93> and using 
the notation of SectionOthe integral becomes 

^ dan n^^,r,x,y) 

The integration path can be taken to run in the upper half plane. Near the origin there 
is a pole, corresponding to the smallest eigenvalue of DkHk- To circle the spectrum 
the integration path would need to pass above the pole but we integrate below the pole 
instead and add the residue K'^°^'^, which we now estimate. 
We see from ( I108> that the zero of dct Vl is at 

p{k) :=ifc- — 
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Figure 3: The integration path for K" 



We have dot n = {t - p(fc))(2ifc + 0{t- p{k)) +0{P)) and 

lim '^^rf^^ -rr^^-^ + Oik). (116) 
p{k) corresponds to Co(fc) — j^'^ + 0{k'^). The terms of Theorem|5lproduce /^p°''= — 

/iP°''= = (e-^'^'* + 0(fce-^''^'*));7+(x)(e-('^+°('^'»ly| - 2kZ+{\y\)), (117) 
< Ce-3'^''*(|fc2+"e-3'=l"-^l| (118) 
+ |/fc2e-^l--y|| + |fc3e-i'=l-l-^l«l|). 

where the C'(-) term comes from (^^} and replacing Co(fc) by l/c'^ (the number ^ is 
just something between j and i). 

To estimate A'''''*', the integral around the rest of the spectrum, we again change to 
T of ( I93> . This leads to integrating above the real axis but below the pole at « ifc as 
in Figure |3l We stay far enough from the pole so that ^ ^^-i remains bounded. This 
allows us to write for small r 



det fc2 + r2 

c.f. (IT08l . Denote R:=%: (C(fc, t") - DkHk)^^. Combining with TheoremlHwe split 
R as follows: 

R — Roo + Rqi + -Rio + -Rii + Ri2 
Roo{x,y) - ^^U+{x){e'^M +2iTZ+{\y\)), 

Roi{x,y)=u;{k,r)U+{x){e'^\y\+2irZ+{\y\)), 

'i(ei^l^+!/l - when xy > 0, 

when xy < 0, 

Ruix,y) = -2TSgii(x(x-y))e-I^^IZ+(|y|) 
|(-a2 + A:2)"iZi2(a;,2/)| < c(|A2"+ie-^MI--«l | 

+ |Ae"5l^-!^l| + |A2e^'^l^l^3lyl|) 

where cj is bounded, /i(T) := minjcImT, 1} for some c > and n G {0, 1}. We 
derived this for small r but by Theorem |6l it actually holds also for large t on our 
integration path. Between small and large r we have Theorems |3] E] and Lemma [Tol 
Everything is continuous on a compact interval, hence bounded. By choosing c appro- 
priately — /i(T) can be used to estimate Re fi2 from above and thus our estimate holds 
everywhere along the integration path. 
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As the explicit terms are analytic we can simply integrate them along the real axis 
(g-C(fc,T")* makes everything small as |t| — > oo). Denote t' := (1 + 2k^)^t. 

= -f7+(x)e-(^-^+'^^)* 1^ p^e-^*' (e'-l^^l + 2iri+(|y|)) ^ 

(119) 



+ (^4fcz(fcVF)-^=^jZ+(|y|) 



where 



1 /"^ 2 

z(a;) := ^ / e"*" dr and /(a;, y) := e^^^ z(x- + y). 

The coefficient of the was obtained from the requirement that K'qq'^{x, y) be con- 
tinuously differentiable in y (because i?oo is)- 



277 

e^^-JF^ - g-^-iF^ ) (120) 



when xy > and zero otherwise. 

/oo 1 

The rest can be estimated by integrating so that Im t > ^k. 

{-dl + k^YK'i^l' = e-^'''*U+{x)0{^^e-i''\y\ + ^e"!^!) (122) 



-dl + fc2)»/i:jf I < Ce-3fe't ( _L_e-f (123) 



\e-h\'^-y\ + l_e-ik\x\-\\y\ 
t 



When k < l/V^ (in addition to fc < e) a slightly different estimate is useful; 
instead of treating the pole separately we integrate around the whole spectrum at once. 
Everything goes essentially as in (m) except that now the integration path goes above 
ifc. This affects the Rqq integral; 



^00 



g(3.^+4fc«)t^+(^) (^/(fcV^^, ^) + f{-kV¥, ^) (124) 
+ (Akz{kVt^) -2k- ^^-L^^Z+{\y\)\ 
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and the estimates for the remainders: 

i-d'^ + k^^Rt = U+ix)Oi-^e-'^ + ie-l^l), (125) 

t ~^ 2 Z 

K-a^ + k^rxill < cl^-^e-"^ (126) 

t ti 

where C depends on c. The other two terms are not effected as they are analytic in r. 



10 Proof of Lemmas [T] and HI 

The above calculations contain more than enough to prove Lemmas ^ and |2] Take 
equal to the e of Section |8l Let us first treat Ka in the case k < Then define 

Ko :== A'^il} + K-^\. (CM and ([125} yield (EoJ and 

= {7+(x)(l + 0(i=(l + M))), 

whence it follows for h £ X that 

£^ e-'^^-^Kfhdy = U+{x)[j e-*-y/i dy + 

= Ut{x) (1(1 + 0(k . y))h dy + Oi}^)) 
= Ut{x)[jhAy + (D{^^)). 
Hence we have estabUshed Mil for k <\/ \/t. 



When l/^/t < k < ko define ATq := A'P°'" + K'^'^ and ATJf^' := A'^^q^' + ATJf . 
Using (II 17t . dl 19> and (I122t from the previous section we get ( I32> and 

Kf' = ?7+(x)(e"^''^'* + + |y|)e-^''-'*), (127) 

A-r = Ut{x)0{^^[l + |j/|)e-^'='*). (128) 
Similarly to the previous case we get 



^rest 
^10 



and \21\ follows. 

Let us then proceed with Ki. Starting again with k < \/^/t write Ki ;= A'" 
Kfl^ + K'fl. ATfo", given by ( I120> . is essentially the first term of (|26j; the difference 
can be absorbed into K2 together with the other two terms, for which see ( I12U and 
( I123> . Thus \2%l is satisfied. ( 13 U is also satisfied. When l/\/t < k < k^ define 
Ki A'P°'" + K[f + Klf + Kfl Using ( fTTSi . ( fT20l . (fT2ll and ivM it is easy to 
see that ( I28> and ( I33> are satisfied. 

Finally, Lemma|3](b) follows from TheoremQand Lemma|3l(a) follows from the 
following theorem: 



33 



Theorem 8. For any e > there exist C, c > such that for any k > e, t > 1 and 

n e {0,1} 

\{e-''^^"''Dl){x,y)\ < Ce~5fc'*-'=l^-!/l 

Proof. Theorem0takes care of things for very large k, say k > r. For k G [e, r] we 
use the Dunford-Cauchy integral 

and estimate the resolvent. The spectrum has k"^ as a lower bound by Lemma|8] Thus 
r can be chosen so that RcC > ■^k* and Rc /ii < Rc /i2 < — c < for some c. 
Asymptotically the path can be chosen to be s ifc^ + se^'" for some a < tt/2. 
The resolvent kernel is then bounded by Ce"'^'^^^': for large \(\ this follows from 
Theorem|6j otherwise just from J104t and continuity. □ 
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